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In [Braz. J. Phys. 30, 27 (2000)] Dickman et al. suggested that self-organized criticality 
can be produced by coupling the activity of an absorbing state model to a dissipation mech- 
anism and adding an external drive. We analyzed the proposed mechanism in [Phys. Rev. E 
73, 025106(R) (2006)] and found that if this mechanism is at work, the finite-size scaling 
found in self-organized criticality will depend on the details of the implementation of dissipa- 
tion and driving. In the preceding comment [Phys. Rev. E XX, XXXX (2008)], Alava et al. 
show that one avalanche exponent in the AS approach becomes independent of dissipation 
and driving. In our reply we clarify their findings and put them in the context of the original 
article. 



PACS numbers: 05.65. -Kb, 05.50.-fq 05.70.Jk, 64.60. Ht 



2 



In [1] we discussed the implications of the absorbing state (AS) mechanism [2] if it were to 
underHe self-organized criticality (SOC). One of the key ingredients of the AS-mechanism is that 
dissipation and driving are made to vanish in the thermodynamic limit. We showed that criticality 
would be reached in the thermodynamic limit for almost any choice of the scaling with system size 
of dissipation and driving (with the effective temperature vanishing, contrary to what is stated in 
[3]). While this choice is thus not important to answer the question if the critical point will be 
reached, it is important when addressing the question how it will be reached. In particular, we 
showed that the critical point would be reached in the limit of slow drive uj — k,> f3/u ("too fast" 
in [1]), but that the observed finite-size scaling exponents would depend on a; — k, and the relative 
correlation length ^/L would vanish asymptotically. 

The present discussion can be phrased in terms of two statements. 1) SOC is universal. 2) The 
AS-mechanism is solely responsible for SOC. In [1] we showed that these are mutually exclusive. 
If the AS-mcchanism is the reason for SOC, then SOC cannot be universal. This would weaken 
SOC considerably since studying models as simple as sandpiles is only sensible if these systems 
show universal behavior. The other possibility is that SOC is universal, which would weaken 
the AS-mechanism, because it does not predict universality. Our analysis was restricted to the 
finite-size scaling of AS observables, such as the order parameter, the correlation length, and the 
susceptibility, whereas the preceding comment [3] by Alava et al. refers to avalanche characteristics. 
In the following we assume that what we found out about the universality of AS observables also 
applies to avalanches. We emphasize that we do not know whether this is true; this assumption is 
made in order to be able to reply to the comment, which makes the same assumption. 

In the following we distinguish between the AS approach intended to explain SOC, and SOC 
itself. Alava et al. do not make this distinction explicit. In the AS approach driving and dissi- 
pation rates are tuned (bulk dissipation as L"'^ and driving as L~'^), whereas in SOC (boundary 
dissipation, driving on a separate time scale) they are set implicitly by the dynamics of the system. 
If the AS approach applies to SOC, then SOC behaviour is obtained within the AS approach by 
taking the thermodynamic limit. Both, our original article [1] as well as the preceding comment 
[3], arc concerned solely with the characterisation of the AS approach. 

It is important to stress that in [1] we did not claim that any exponents, neither those describing 
avalanches nor those characterising the activity, in standard SOC models are non-universal. This 
deserves clarification, because the opposite is stated in the abstract of [3] . As discussed in [1] , there 
are instances of SOC exponents being identical across a wide range of models [4, 5, 6], while others 
are not [5, 6]. 
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We do claim (with the proviso stated above), however, that avalanche size exponents would be 
non-universal if the AS mechanism [2] was applicable to SOC [1]. Alava et al. challenge this claim 
by showing that within the AS approach Tg is independent of the scaling of driving and dissipation 
in the slow driving limit. This limit corresponds to the separation of time scales in SOC. 

In [1] we explicitly mention avalanche exponents only once: "In the AS approach also the 
avalanche size exponents show a clear, immediate dependence on the choice of the two exponents k 
and a;." It is correct that the avalanche size exponents Tg and Dl (see below) depend on k and a;, 
but as Alava et al. show, for a; large enough Tg becomes independent of k and a;. In this case, as 
we show below, nonetheless depends on n. The non-universality of within the AS approach 
implies that it does not explain universal SOC. This is the same conclusion we reached by studying 
AS observables (order parameter, correlation length, susceptibility). 

We agree that the derivation in [1] necessarily breaks down in the large u limit. We explicitly 
assumed finite bounds for both, k and cj, and we did not discuss the case of uj being greater than 
the dynamical exponent, which is the regime studied by Alava et al.. We agree that this regime is 
the most important one for SOC. 

Alava et al. have chosen an observable that is independent of the external drive for sufficiently 
large a;, but its finite-size scaling turns out to depend on the scaling of the dissipation. This can be 
seen by a finite-size scaling analysis of the characteristic avalanche size Sc- According to Ref. [3], 
Eq. (3), Sc Of , where is some "cut-off scale". Usually, the exponent D is reserved for the finite 
size scaling of Sc, also known as the avalanche dimension [5], which we call Dl in the following, 
so that Sc oc L^^. Combining Eq. (4) and Eq. (3) of [3], they find Dl = n/{2 — Tg) with Tg being 
independent of external drive and dissipation. This is a very surprising rcsTilt, because regardless 
of whether or not the AS-mechanism applies, the avalanche dimension Dl is deeply rooted in the 
model and in general directly related to the field theory of the corresponding depinning transition 
[5, 7, 8, 9]. There are several models [4, 5, 6] which display a universal avalanche dimension Dl 
and an avalanche size exponent Tg = 2 — ^\/Dl which depends on the details of the driving of the 
model [5], with the first moment scaling like (s) oc LPi'^ . For these conservative models, it is straight 
forward to devise a method to produce any exponent Tg in the interval [2 — 2/1)^,2 — 1/Dl] by 
effectively tuning 71. This can be achieved by changing the driving mechanism, [10] as the driving 
mechanism leaves Dl unchanged. 

We would have expected that changing 71 by introducing dissipation would have the same 
effect, i.e. varying Tg and constant Dl. However, Alava et al. find that Dl depends on k, while 
Tg remains unchanged, a very interesting numerical finding we do not dispute. It has thus been 
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established that, under the appropriate conditions, both Tg and Dl can be tuned. Importantly, 
Dl can be tuned in the SOC-regime of the AS approach (large uj) by changing k. 

Inasmuch as avalanche exponents pertain to the discussion the universality of Tg within the AS 
approach supports the case for universal SOC being generated by the AS mechanism. But, the 
implicit finding by Alava et al. that Dl depends on n confirms that, apparently, the AS approach 
does not produce universal SOC. 

In order to address the question whether the AS-mechanism is operating in SOC models (uni- 
versal or not), one needs to probe its presence either directly or test its implications. In [1] we 
have shown that the AS-mechanism would (almost always) imply a vanishing relative correlation 
length ^/L and a finite-size scaling of the AS order parameter, characterized by exponents and 
7//X, that would depend on the scaling of dissipation and drive, parameterized by k and u. We 
stated explicitly how f3/fi and ^ / depend on k and a;, while further analysis is necessary for the 
dependence of the avalanche exponents on k and uj. At the present stage, a more promising route 
than studying avalanches therefore seems to be the study of AS observables in SOC systems. 
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